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Coherent-scatter computed tomography 



5 

The present invention relates to the field of coherent-scatter computed 
tomography (CSCT), where a fan-beam is applied to an object of interest In particular, 
the present invention relates to a data processing device for performing a reconstruction 
of computed tomography data, to a computed tomography apparatus for examination of 
10 an object of interest, to a method of performing a reconstruction of computed 

tomography data and to a computer program for a data processor for performing a 
reconstruction of computed tomography data. 



15 US 4,751,722 describes a device based on the principle of registration of 

an angled distribution of coherent scattered radiation within angles of 1° to 12° as 
related to the direction of the beam. As set forth in the US 4,751,722, the main fraction 
of elastic scattered radiation is concentrated within angles of less than 12°, and the 
scattered radiation has a characteristic angle dependency with well marked maxima, the 
positions of which are determined by the irradiated substance itself. As the distribution 
of the intensity of the coherently scattered radiation in small angles depends on 
molecular structure of the substance, different substances having equal absorption 
capacity (which cannot be differentiated with conventional transillumination or CT) can 
be distinguished according to the distribution of the intensity of the angled scattering of 
coherent radiation typical for each substance. 

Due to the improved capabilities of such systems to distinguish different 
object materials, such systems find more and more application in medical or industrial 
fields. 

The dominant component of low-angle scatter is coherent scatter. 
Because coherent scatter exhibits interference effects which depend on the atomic 
arrangement of the scattering sample, coherent-scatter computed tomography (CSCT) is 
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in principle a sensitive technique for imaging spatial variations in the molecular 
structure of tissues across a 2D object section. 

Harding et al. **Energy-dispersive x-ray diffraction tomography" Phys. 
Med. Biol., 1990, Vol. 35, No. 1, 33-41 describes a tomographic imaging technique 
5 based on an energy analysis at fixed angle, of coherent x-ray scatter excited in an object 
by polychromatic radiation. According to this method, a radiation beam is created by the 
use of suitable aperture systems, which has the form of a pencil and thus is also referred 
to as a pencil beam. Opposite to the pencil beam source, one detector element suitable 
for an energy analysis is arranged for detecting the pencil beam altered by the object of 
10 interest. 

A coherent scatter set-up applying a fan-beam primary beam and a 2D 
detector in combination with CT was described in US 6,470,067 Bl. The shortcoming 
of the angle-dispersive set-up in combination with a polychromatic source are blurred 
scatter functions, which is described in e.g. Schneider et al. "Coherent Scatter 
15 Computed Tomography Applying a Fan-Beam Geometry" Proc. SPIE, 2001, Vol. 4320 
754-763. 

To become a competitive modality in the fields of medical imaging or 
non-destructive testing, the implemented reconstruction algorithm should assure both 
good image quality and short reconstruction times. 

So far, the projection data acquired with fan-beam CSCT is reconstructed 
with the help for example, algebraic reconstruction techniques (ART), since ART has 
been shown to be highly versatile, for example, by J. A. Grant et al. "A reconstruction 
strategy suited to x-ray diffraction tomography", J.Opt. Soc. Am A12, 291-300 (1995). 
However, due to the computational complexity of such iterative reconstruction, such 
methods require a relatively long reconstruction time. 

Furthermore, all known approaches apply a non-exact reconstruction 

technique. 

It is an object of the present invention to provide for an improved 
reconstruction of tomography data. 

According to exemplary embodiment of the present invention as set forth 
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in claim 1, the above object may be solved by a data processing device for performing a 
reconstruction of computed tomography data, wherein the computed tomography data is 
reconstructed from acquired CT (computed tomography) data, comprising at least a 
partial spectrum acquired by using a detector comprising energy resolving detector 
5 elements. Furthermore, a memory is provided for storing at least one of the acquired CT 
data and the computed tomography data. Also, a processor is provided, which is adapted 
to determine a wave-vector transfer by using the at least partial spectrum and to 
determine a reconstruction volume. A dimension of the reconstruction volume is 
determined by the wave-vector transfer. The wave-vector transfer represents curved 

10 lines in the reconstruction volume. Then, according to an aspect of the present 

invention, a rearranging of the CT data is performed, such that it corresponds to an 
acquisition along a desired source trajectory in the reconstruction volume. 

Advantageously, the data processing device according to this exemplary 
embodiment of the present invention may perform a quasi-exact reconstruction 

15 technique for reconstruction data, such as an image, from read-outs of the detector. In 
this context, the term "reconstruction" is to be understood as covering steps or processes 
from the reading of the measured data from the detector to the display of the actually 
reconstructed data on, for example, a display or an output of information such as a 
material discrimination based on these read-outs. 

20 The interpretation of the projection data as curved lines or line integrals 

in the reconstruction space and the appropriate resorting of the acquired CT data such 
that it corresponds to an acquisition, where the source of radiation used is displaced 
along a desired source trajectory, may allow for realizing exact helical reconstruction 
algorithms. 

25 In other words, according to this exemplary embodiment of the present 

invention, the source of radiation may be displaced along a first trajectory, such as a 
circle, for data acquisition. Then, the data is rearranged such that it corresponds to data 
as if it had been acquired with a source of radiation displaced along a second trajectory 
which may be different to the first trajectory. The second trajectory may be a helix. 

30 According to another exemplary embodiment of the present invention as 

set forth in claim 2, the acquired CT data is acquired with a first source movement and 
the data is rearranged in the reconstruction volume such that it corresponds to data 
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acquired with a second, different source movement. E.g. the first source movement may 
be a circle and the second source movement may be a helix. Due to this reordering or 
reorganization of the acquired CT data, the exact helical reconstruction algorithms, 
which are known in the art as described e.g. in Katsevich "Analysis of an exact 
5 inversion algorithm for spiral cone-beam CT", Phys. Med. Biol., vol. 47, p 2583-2597, 
2002, which is hereby incorporated by reference, may be used for quasi-exactly 
reconstructing the data. 

According to another exemplary embodiment of the present invention as 
set forth in claim 3, a filtered back-projection is performed along the curved lines, 

1 0 which may be hyperbolas in the reconstruction volume. 

According to another exemplary embodiment of the present invention, 
the three dimensions of the reconstruction volume are defined by the wave-vector 
transfer and by two linear independent vectors of the rotation plane of the source of 
radiation. Furthermore, the detector is a two-dimensional detector, which allows to 

15 achieve the full energy spectrum of the source of radiation, for example, of a 

polychromatic x-ray source. This yields a relatively wide spectrum of the wave-vector 
transfers of the scattered x-ray photons. Due to this, redundant data is acquired, which ■ 
may be used to obtain a better image quality. 

According to another exemplary embodiment of the present invention as 

20 set forth in claim 5, the rearranging of the acquired CT data is performed by using 

John's equation. According to another exemplary embodiment of the present invention 
as set forth in claim 6, a very simple linear equation is provided, which fulfills the data 
acquisition of a helical trajectory in the reconstruction volume, allowing for a fast and 
accurate and quasi-exact reconstruction of the computed tomography data. 

25 According to another exemplary embodiment of the present invention as 

set forth in claim 7, a computed tomography apparatus is provided, with a scatter 
radiation detector arranged at a detector unit opposite to an x-ray source, with an offset 
with respect to a slice plane of a fan-shaped x-ray beam generated by the x-ray source in 
a direction parallel to a rotational axis of the x-ray source and the scatter radiation 

30 detector. The scatter radiation detector comprises a plurality of energy resolving 

detector elements. According to an aspect of this exemplary embodiment of the present 
invention, a rearranging of the acquired CT data acquired by using the scatter radiation 



WO 2005/036467 



5 



PCT/IB2004/051968 



detector is performed such that it corresponds to an acquisition, where the x-ray source 
is displaced along a desired source trajectory in the reconstruction volume. In other 
words, the acquired CT data is rearranged as if it had been acquired by means of a 
desired movement of the scatter radiation detector and the x-ray source which may be 
5 different to the actual movement of the x-ray source during the data acquisition. 

Advantageously, according to this exemplary embodiment of the present 
invention, a computed tomography apparatus is provided, which may perform a quasi- 
exact helical reconstruction of, for example, an x-ray image. Advantageously, this may 
allow for improved image quality and for a fast reconstruction. 

According to another exemplary embodiment of the present invention as 
set forth in claim 8, the scatter radiation detector is a two-dimensional detector, i.e. a 2D 
energy resolving detector, which may allow to achieve the full energy spectrum of, for 
example, a polychromatic x-ray source. This may yield a relatively wide spectrum of the 
wave-vector transfers of the scattered x-ray photons. An interpretation of the projection 
data as line integrals in the reconstruction volume and the resorting of the data as if it had 
been taken along the helical source trajectory may allow for applying exact helical 
reconstruction algorithms. Furthermore, redundant data may be used to obtain a better 
image quality. 

Claims 9 to 1 1 provide for further exemplary embodiments of the 
computed tomography apparatus according to the present invention. 

According to another exemplary embodiment of the present invention as 
set forth in claim 12, a method of performing a reconstruction of computed tomography 
data is provided, comprising a rearranging of acquired CT data, acquired by means of an 
energy resolving detector, which may be a two-dimensional energy resolving detector, 
as if it had been taken by an acquisition along a helical source trajectory. 

Advantageously, according to this exemplary embodiment of the present 
invention, a very fast and efficient method may be provided, allowing for an exact 
reconstruction of the data, which may allow for an improved image quality. 

Claims 13 and 14 provide for further exemplary embodiments of the 
method according to the present invention. 

According to another exemplary embodiment of the present invention as 
set forth in claim 1 5, a computer program for a data processor for performing a 
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reconstruction of computed tomography data is provided. The computer program 
according to the present invention is preferably loaded into a working memory of the 
data processor. The data processor is thus equipped to carry out the method of the 
invention. The computer program may be stored on a computer readable medium, such 
5 as a CD-ROM. The computer program may also be presented over a network such as 
the Worldwide Web, and can be downloaded into the working memory of a data 
processor from such a network. 

It may be seen as the gist of an exemplary embodiment of the present 
invention that a resorting or extrapolation of the data acquired during a circular 
10 acquisition (i.e. the data read-outs from, for example, a two-dimensional energy 

resolving detector) is performed, such that the data corresponds to an acquisition along 
a different source trajectory as actually used for acquisition in the x-y-q space, or the 
reconstruction volume. In other words, e.g. data acquired during a circular acquisition is 
resorted and/or extrapolated such that it represents data in an x-y-q space, as if it had 
1 5 been taken by an acquisition with a helical source trajectory. Advantageously, this may 
allow for a quasi-exact helical reconstruction, since known reconstruction algorithms 
adapted to helical source trajectories may be used. Furthermore, according to an aspect 
of the present invention, a two-dimensional energy resolving detector may be used, 
allowing to achieve a full energy spectrum of the source of radiation, which yields a 
relatively wide spectrum of the wave-vector transfers of the scattered photons. 
Advantageously, this may allow for an improved image quality and for an improved 
data reconstruction. 

These and other aspects of the present invention will become apparent 
from and elucidated with reference to the embodiments described hereinafter. 

Exemplary embodiments of the present invention will be described in the 
following, with reference to the following drawings: 

Fig. 1 shows a schematic representation of an exemplary embodiment of 
a CSCT scanner according to the present invention. 
Fig. 2 shows a schematic representation of the geometry of the CSCT 
scanner of Fig. 1 for the measurement of coherent scatter radiation. 
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Fig. 3 shows another schematic representation of the geometry of the 
CSCT scanner of Fig. 1. 

Fig. 4 shows another schematic representation of the measurement 
geometry of the CSCT scanner of Fig. 1 for further explaining the present 
5 invention. 

Fig. 5 shows a schematic representation of a side view of the geometry of 
the CSCT scanner of Fig. 1. 

Fig. 6 shows a schematic representation of a multi-line CSCT scanner 
according to another exemplary embodiment of the present invention. 
10 Figs. 7a and 7b are schematic drawings to explain the idea of John's 

equation. 

Fig. 8 shows characteristics of the wave-vector transfer q in dependence 
on the distance D between the scatter center and the detector for further 
explaining the present invention. 
15 Fig. 9 shows a simplified schematic representation of an exemplary 

embodiment of a data processing device according to the present 
invention. 



20 Fig. 1 shows an exemplary embodiment of a computer tomograph 

according to the present invention. With reference to this exemplary embodiment, the 
present invention will be described for the application in baggage inspection to detect 
hazardous materials such as explosives in items of baggage. However, it has to be noted 
that the present invention is not limited to applications in the field of baggage 

25 inspection, but can also be used in other industrial or medical applications, such as for 
example in bone imaging or a discrimination of tissue types in medical applications. 

The apparatus depicted in Fig. 1 is a fan-beam CSCT scanner, which 
allows in combination with an energy-resolving detector and with tomographic 
reconstruction a good spectral resolution, even with a polychromatic primary fan-beam. 

30 The CSCT scanner depicted in Fig. 1 comprises a gantry 1, which is rotatable around a 
rotational axis 2. The gantry 1 is driven by means of a motor 3. Reference character 4 
designates a source of radiation, such as an x-ray source, which, according to and aspect 



WO 2005/036467 



8 



PCT7IB2004/051968 



of the present invention, emits a polychromatic radiation. 

Reference character 5 designates a first aperture system, which forms the 
radiation beam emitted from the radiation source 4 to a cone shaped radiation beam 6. 
Furthermore, there is provided another aperture system 9 consisting of a diaphragm or a 
5 slit collimator. The aperture system 9 has the form of a slit 10, such that the radiation 
emitted from the source of radiation 4 is formed into a fan-beam 1 1 . According to a 
variant of this exemplary embodiment of the present invention, the first aperture system 
5 may also be omitted and only the second aperture system 9 may be provided. 

The fan-beam 1 1 is directed such that it penetrates the item of baggage 7, 

10 arranged in the center of the gantry 1, i.e. in an examination region of the CSCT scanner 
and impinges onto the detector 8. As may be taken from Fig. 1, the detector 8 is 
arranged on the gantry 1 opposite to the radiation source 4, such that the slice plane of 
the fan-beam 1 1 intersects a row or line 15 of the detector 8. The detector 8 depicted in 
Fig. 1 has seven detector lines, each comprising a plurality of detector elements. As 

15 mentioned above, the detector 8 is arranged such that the primary radiation detector 15, 
i.e. the middle line of the detector 8 is in the slice plane of the fan-beam 11. 

As can be taken from Fig. 1, the detector 8 comprises two types of 
radiation detector lines: a first type of detector lines 30 and 34, which are indicated 
without hatching in Fig. 1, which are detector lines consisting of energy resolving 

20 detector cells. According to an aspect of the present invention, these first detector 
elements (lines 30 and 34) are energy-resolving detector elements. Preferably, the 
energy resolving detector elements are direct-converting semiconductor detectors. 
Direct-converting semiconductor detectors directly convert the radiation into electrical 
charges - without scintillation. Preferably, these direct-converting semiconductor 

25 detectors have an energy resolution better than 10 % FWHM, i.e. AE/E< 0.1, with AE 
being the full-width at half maximum (FWHM) of the energy resolution of the detector. 

Such detector cells of lines 30 and 34 may be cadmiumtelluride or 
CdZnTe (CZT) based detector cells, which are both outside of the slice plane of the fan- 
beam 1 1. In other words, both energy resolving lines 30 and 34 are arranged at the 

30 gantry 1 opposite to the x-ray source 4 with an offset from the slice plane in a direction 
parallel to the rotational axis 2. The detector line 30 is arranged with a positive offset 
with respect to the direction of the rotational axis 2 depicted in Fig. 1, whereas the line 
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34 is arranged with a negative offset from the slice plane with respect to the direction of 
the rotational axis 2 depicted in Fig. 1. 

The detector lines 30 and 34 are arranged at the gantry 1 such that they 
are parallel to the slice plane and out of the slice plane with an offset in a positive or 
5 negative direction of the rotational axis 2 of the gantry 1, such that they receive or 

measure a scatter radiation scattered from the item of baggage 7 in the examination area 
of the CSCT scanner. Thus, in the following, lines 30 and 34 will also be referred to as 
scatter radiation detector. It has to be noted that according to a preferred embodiment of 
the present invention, the energy resolving lines 30 and 34 form a two dimensional 

10 energy resolving scatter radiation detector. In other words, preferably, a plurality of 
energy resolving lines 30 and 34 are arranged on both sides of the line 15. Thus, a 
detector 8, having only a non energy resolving line 15 in the slice plane and having 
energy resolving lines arranged spaced from the slice plane, may be advantageous. 
Thus, if, in the following the term "scatter radiation detector" is used, it includes any 

1 5 detector with at least two lines of energy resolving detector cells, which are arranged out 
of the fan plane of the fan-beam 11, such that it receives photons scattered from the item 
of baggage 7. 

The second type of detector lines provided on the detector 8, which are 
indicated by a hatching, are scintillator cells. In particular, line 15 is arranged such that 

20 it is in the slice plane of the fan-beam 1 1 and measures the attenuation of the radiation 
emitted by the source of radiation 4, caused by the item of baggage 7 in the examination 
area. As depicted in Fig. 1, right and left of the line 15, there may be provided further 
detector lines including scintillator detector cells. However, as indicated above, there 
may be only one detector line with scintillator cells, namely the line 15 in the slice 

25 plane, and the remaining lines of the detector 8 may be energy resolving lines. By this, a 
two dimensional energy resolving detector (such a scatter radiation detector) may be 
provided. 

As already indicated with respect to the energy resolving lines 30 and 34, 
the provision of only the line 15 measuring the attenuation caused by the item of 
30 baggage 7 of the primary beam of the fan-beam 1 1 in the slice plane may be sufficient 
However, as in the case of the energy resolving lines 30 and 34, a provision of a 
plurality of detector lines 32, each comprising a plurality of scintillator cells, may 



WO 2005/036467 



10 



PCT/IB2004/051968 



further increase the measurement speed of the CSCT scanner. In the following, the term 
"primary radiation detector" will be used to refer to a detector, including at least one 
line of scintillator or similar detector cells for measuring an attenuation of the primary 
radiation of the fan-beam 1 1 . 
5 According to a preferred embodiment of the present invention, a plurality 

of energy resolving lines 30 and 34 are provided, such that a two-dimensional energy 
resolving detector is provided. This may allow to gather additional, redundant data, 
which may be used later on for improving the image quality. Furthermore, this may 
allow to reduce a scanning time. 

10 With one energy-resolving line, only a fraction of the full spectrum for 

wave-vector transfers may be obtained. In addition, the spectrum may have a lower 
boundary, caused by a distance of the detector line to the slice plane or scan plane. 
Thus, using a 2D (two dimensional) detector may lead to a wider spectrum of the wave- 
vector transfers. In particular, the lower boundary may be shifted to relatively small 

1 5 values. Furthermore, a 2D detector may allow to yield redundant data. By pre- 
processing this data in accordance with an aspect of an exemplary embodiment of the 
present invention, it can be used as input for an exact reconstruction technique. 

As may be taken from Fig. 1, the detector cells of the detector 8 are 
arranged in lines and columns, wherein the columns are parallel to the rotational axis 2, 

20 whereas the lines are arranged in planes perpendicular to the rotational axis 2 and 
parallel to the slice plane of the fan-beam 11. 

The apertures of the aperture systems 5 and 9 are adapted to the 
dimensions of the detector 8 such that the scanned area of the item of baggage 7 is 
within the fan-beam 1 1 and that the detector 8 covers the complete scanning area. 

25 Advantageously, this allows to avoid unnecessary excess radiation applied to the item of 
baggage 7. During a scan of the item of baggage 7, the radiation source 4, the aperture 
systems 5 and 9 and the detector 8 are rotated along the gantry 1 in the direction 
indicated with arrow 16. For rotation of the gantry 1 with the source of radiation 4, the 
aperture systems 5 and 9 and the detector 15, the motor 3 is connected to a motor 

30 control unit 17, which is connected to a calculation unit 18. 

In Fig. 1, the item of baggage 7 is disposed on a conveyor belt 19. 
According to an exemplary embodiment of the present invention, as mentioned above, a 
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circular data acquisition is performed where the x-ray source 4 is displaced along a 
circular source trajectory, i.e. is rotated in a rotational plane around the rotation axis 
without a movement of the item of baggage 7. Thus, the conveyor belt 19 is operated 
such that the item of baggage is immobile when the x-ray source 4 performs a circular 
5 movement (i.e. a circular acquisition) around the item of baggage 7 allowing for a 
sufficient data acquisition. Then, the item of baggage 7 is moved by e.g. a preset 
increment, and then another scan of another slice of the item of baggage 7 is performed 
during which the conveyor belt 19 is not operated. 

The detector 8 is connected to a calculation unit 18. The calculation unit 

10 18 receives the detection results, i.e. the readouts from the detector elements of the 

detector 8 and determines a scanning result on the basis of the scanning results from the 
detector 8, i.e. from the energy resolving lines 30 and 34 and the lines 15 and 32 for 
measuring the attenuation of the primary radiation of the fan-beam 1 1 . In addition to 
that, the calculation unit 18 communicates with the motor control unit 17 in order to 

15 coordinate the movement of the gantry 1 with the motors 3 and 20 or with the conveyor 
belt 19. 

The calculation unit 18 is adapted for reconstructing an image from 
readouts of the primary radiation detector, i.e. detector lines 15 and 32 and the scatter 
radiation detector, i.e. lines 30 and 34. The image generated by the calculation unit 18 

20 may be output to a display (not shown in Fig. 1) via an interface 22. 

The calculation unit, which may be realized by a data processor, may be 
adapted to perform a filtered back-projection on the read-outs from the detector element 
of the detector 8, i.e. from the read-outs from the energy resolving lines 30 and 34 and 
the lines 15 and 32 for measuring the attenuation of the primary radiation of the fan- 

25 beam 1 1 . The back-projection performed in the calculation unit 1 8, which forms part of 
the image reconstruction will be described in further detail with reference to Fig. 7. 

Furthermore, the calculation unit 18 may be adapted for the detection of 
explosives in the item of baggage 7 on the basis of the readouts of the lines 30 and 34 
and 15 and 32. This can be made automatically by reconstructing scatter functions from 

30 the readouts of these detector lines and comparing them to tables including 
characteristic measurement values of explosives determined during preceding 
measurements. In case the calculation unit 18 determines that the measurement values 
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read out from the detector 8 match with characteristic measurement values of an 
explosive, the calculation unit 18 automatically outputs an alarm via a loudspeaker 21 . 

In particular, the calculation unit 18 may by adapted to perform a 
determination of the wave-vector transfer by using at least a partial spectrum acquired 
5 from the read-outs from the detector 8 and, in particular, from the read-outs of the 
energy resolving lines 30 and 34. As may be taken from Fig. 1, the data is acquired 
during a circular acquisition, referring to the source trajectory in the x-y-z space 

Then, the measured data is interpreted as line integrals in the x-y-q space, 
where q is the wave-vector transfer and x and y are linear independent vectors in the 

10 rotation plane of the source of radiation 4. Then, the measured data is resorted and 
extrapolated, such that it corresponds to an acquisition along a helical trajectory in the 
x -y-q space, i.e. it is reordered to be as if it had been taken along a helical source 
trajectory. Then, the data may be pre-processed in order to allow that the data is applied 
to conventional helical reconstruction algorithms, for example, the reconstruction 

1 5 technique described by Katsevich "Analysis of an exact inversion algorithm for spiral 
cone-beam CT", Phys. Med. Biol., vol. 47, p. 2583-2597, 2002, which is hereby 
incorporated by reference. Finally, the data is back-projected. This back-projection may 
be performed along the curved lines, such as hyperbolas in the x-y-q space of the wave- 
vector transfer. This operation will be described in further detail with reference to Fig. 

20 6. 

During the subsequent description of Figs. 2 to 5, the same reference 
numbers as used in Fig. 1 will be used for the same or corresponding elements. 

Fig. 2 shows a simplified schematic representation of a geometry of the 
CSCT scanning system depicted in Fig. 1. As may be taken from Fig. 2, the x-ray source 
25 4 emits the fan-beam 1 1 such that it includes the item of baggage 7 in this case having a 
diameter of u and covers the entire detector 8. The diameter of the object region may, 
for example, be 100 cm. In this case, an angle a of the fan-beam 1 1 may be 80 °. In such 
an arrangement, a distance v from the x-ray source 4 to the center of the object region is 
approximately 80 cm and the distance of the detector 8, i.e. of the individual detector 
30 cells from the x-ray source 4 is approximately w = 150 cm. 

As can be taken from Fig. 2, according to an aspect of the present 
invention, the detector cells or lines can be provided with collimators 40 to avoid that 
the cells or lines measure unwanted radiation having a different scatter angle. The 
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collimators 40 have the form of blades or lamellas, which can be focused towards the 
source. The spacing of the lamellas can be chosen independently from the spacing of the 
detector elements. 

Instead of a bent detector 8 as depicted in Figs. 1 and 2, it is also possible 
5 to use a flat detector array. 

Fig. 3 shows another schematic representation of a detector geometry as 
used in the CSCT scanner of Fig. 1. As already described with reference to Fig. 1, the 
detector 8 may comprise one, two or more energy resolving detector lines 30 and 34 and 
a plurality of lines 15 and 32 for measuring the attenuation of the primary fan-beam 
10 caused by the item of baggage 7. As may be taken from Fig. 3, preferably the detector 8 
is arranged such that one line of the lines 15 and 32, preferably the middle line 15 of the 
detector 8, is within the slice plane of the fan-beam 1 1 and thereby measures the 
attenuation in the primary radiation. As indicated by arrow 42, the x-ray source 4 and 
the detector 8 are rotated together around the item of baggage in the rotation plane to 
1 5 acquire projections from different angles to perform a circular acquisition. 

As depicted in Fig. 3, the detector 8 comprises a plurality of columns L 
Fig. 4 shows another schematic representation of the geometry of the 
CSCT scanner depicted in Fig. 1 for further explaining the present invention. In Fig. 4, a 
detector 46 is depicted, comprising only one line 15 and only one line 30. The line 15 is 
20 arranged in the slice plane of the fan-beam 1 1 formed by the aperture system 9, which in •• 
this case is a slit collimator and generated by means of the source of radiation or x-ray 
source 4. The line 15 comprises, for example, scintillator cells or other suitable cells for 
measuring the attenuation of the primary beam of the fan-beam 1 1 and allows for an 
integral measurement of the attenuation of the primary fan-beam caused by the object of 
25 interest in the object region or examination region. 

Line 30 depicted in Fig. 4 includes energy resolving cells. As may be 
taken from Fig. 4, the line 30 is arranged parallel to the slice plane of the fan-beam 1 1 
but out of the plane. In other words, the line 30 is arranged in a plane parallel to the 
slice plane and parallel to the line 15. 
30 Reference numeral 44 indicates a scatter radiation, Le. a photon scattered 

by the object of interest, such as the item of baggage. As may be taken from Fig. 4, the 
scatter radiation leaves the slice plane and impinges onto a detector cell of the line 30. 
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Fig. 5 shows a side view of the detector geometry of the CSCT scanner 
of Fig. 1. Fig. 5 can also be contemplated as showing a side view of Fig. 4, where, 
however, instead only the provision of one line 30 and one line 15, in Fig. 5, there is 
provided a plurality of detector lines 32 between the line 30 and the line 15. The 
5 detector element D$ of the line 30 is an energy resolving detector element. The detector 
element D% is arranged with a fixed distance a from the slice plane of the primary fan- 
beam. According to an aspect of the present invention, for each detector element D% of 
the column t and for each projection <D (see Fig. 3) an intensity spectrum / (E, t, <D) is 
measured. Performing this measurement for a plurality of projections <D along a circular 

10 scan path, a three-dimensional dataset is acquired. Each object pixel is described by 
three coordinates (x, y, q). Thus, according to an aspect of the present invention, for 
reconstructing an image or for reconstructing further information from the three- 
dimensional dataset, a 3D — > 3D reconstruction method can be used. 

Fig. 6 shows a schematic drawing of an exemplary embodiment of a 

15 multi-line CSCT scanner. This scanner is provided with a detector 48, comprising a 
plurality of lines of energy resolving detector elements, which may be the same as the 
ones described with reference to Fig. 1 . The source of radiation 9 is provided with 
collimator means, such that it generates a fan-beam of x-rays. The arrangement of the 
detector 48 and the source of radiation 49 is such that the detector 48 is focus centered. 

20 The view depicted in Fig. 6 is parallel to the scanned plane or slice plane in order to 
further clarify the scanning process out of the x-y plane, i.e. the rotation plane of the 
source of radiation 49 and the detector 48. As may be taken from Fig. 6, a distance 
between the source of radiation 49 and the detector 48 is indicated as "SD"; a distance 
between the source 49 and the center of rotation 47 is indicated by S, a distance between 

25 the scatter center and the detector 48 is indicated by d, a distance between a detector 
element receiving radiation and the scanned plane or sliced plane is given by a and h 
indicates a height of the detector 48. 

The z coordinate axis is a normal on the center of the rotation plane of 
the source of radiation 49, i.e. the axis of rotation of the source of radiation 49. Thejy 

30 coordinate is in the rotation plane of the source of radiation. 

As may be taken from Fig. 6, for the following description, a CSCT 
scanner with, for example, a polychromatic x-ray source 49 and a D detector 49 is 
considered. The detector comprises or consists of energy resolving detector elements, 
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which may be similar to the ones described with reference to Fig. 1 . The emitted x-rays 
have been collimated, such that a fan-beam radiates the object of interest located in an 
area around the center of radiation 47. 

According to an exemplary embodiment of the present invention, the 
5 following method of operation may be applied in the above scanner, or in the scanner 
described with reference to Fig. 1, for reconstructing the CSCT data, i.e. for 
reconstructing an image from the read-outs of the detectors 8 and 48. 

Step 1: The data is measured during a circular acquisition referring to the 
source trajectory in the x-y-z space. In other words, read-outs are gathered from 
10 detectors 8 or 48, while the source of radiations 4 and 49 and the detectors 8 and 48 are 
rotated around the object of interest in a rotational plane. The read-outs are referred to 
as measured data or acquired CSCT data. The measured CSCT data are interpreted as 
line integrals in the x-y-q space, where q represents the wave-vector transfers. The 
calculation of the wave-vector transfers will be described later on. 
15 Step 2: The acquired CSCT data is resorted and extrapolated, such that it 

corresponds to an acquisition along a helical trajectory in the x-y-q space. 

Step 3: A further step may be performed in order to pre-process the data 
according to conventional helical reconstruction algorithms, such as, for example, the 
exact reconstruction technique described by Katsevich "Analysis of an exact inversion 
20 algorithm for spiral cone-beam CT", Phys. Me<L Biol., vol. 47, p. 2583-2597, 2002, 
which is hereby incorporated by reference. 

Step 4: Then, the resorted and/or extrapolated data may be back 
projected. This back-projection may be performed along the curved lines in the x-y-q 
space. These curved lines may, for example, be hyperbolas. 
25 This operation, in particular step 2, will be described in more detail in the 

following: 

CSCT makes use of coherently scattered x-rays, in order to reconstruct 
the coherent scatter form factor F^ig). The differential cross-section for coherently 
scattered x-rays da^^ /dQ is given by 

where r e denotes the classical electron radius, and 0 the angle between 
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the incoming and the scattered x-rays. The wave-vector transfer q causing the deviation 
of the photon by the angle 0 is defined by 

? = — sin(©/2), (2) 

with the energy E of the corresponding x-ray photon, Planck's constant h 
5 and the speed of light c. For scattering under small angles, e.g. the angle regime of 
interest here is between 0 and 5°, sin(0/2) can be approximated by 0/2, and Eq.(2) 
can be written as 

E © 

»-CT (3) 

According to Fig. 6, the scatter angle is given by the distance d of the 
1 0 scatter center from the detector and the distance a of the detector element which 
receives the scattered radiation from the scanning plane: 

tan©«© = ~. (4) 

Together with Eq. (3), this yields: 
E a 

15 The x-y-q space, Eq. (5) describes hyperbolas. These hyperbolas can be 

approximated by straight lines. Among several possibilities is, for example, an 
approximation such that the area under a straight line matches the area of the 
corresponding hyperbola. Another approximation is described here. The straight line 
intersects the hyperbola at the beginning (d^ ) and at the end {d^ ) of the region of 

20 interest: 



q 2hc 



d mm Xd m*x 



aE. (6) 



A scanning system is considered, such as the ones depicted in Figs. 1 and 
6, where the real path can be resorted and extrapolated such that they correspond to an 
acquisition along a helical trajectory in the x-y-q space. The extrapolation of the data 
25 measured on a circular trajectory to virtual neighboring trajectories is done by using 
John ! s Equation. The idea of John's Equation is that the space of line integrals through 
3D space is 4D, thus the mapping from object function to its line integral function 
creates an additional dimension, as described in S. K. Patch, "Consistency conditions 
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upon 3D CT data and the wave equation", Phys. Med. Biol. 47, 2637-2650. 
US Patent 6173030 (1999), which is hereby incorporated by reference. 

Figs. 7a and 7b are schematic drawings in order to further explain the 
idea of John's Equation. The line integrals in Fig. 7b for the virtual source position q are 
5 extrapolated from the line integrals measured for source positions for q = 0 shown in 
Fig. 7a. 

This implies redundancy in line integral space, which is then used to 
construct unmeasured from measured data as described in S. K. Patch "Computation of 
unmeasured third-generation VCT views from measured views", IEEE Trans. Med. 

10 Img. MI-21, 801-813. US Patent 6292526 (1999), which is hereby incorporated by 
reference. John's Equation is parametrized as follows for our set of geometry 
parameters, as described in M. Defrise, F. Noo, H. Kudo, "Improved 2D rebinning of 
helical cone-beam CT data using John's equation", Proc. 2002 IEEE Nuclear Science 
and Medical Imaging Symposium, Norfolk (VA), Paper M10-74, which is hereby 

1 5 incorporated by reference: 

R 2 gu< ~2ug a - OR 2 + u 2 )g ua -R gaa -uag aa = 0, (7) 

where R is the distance from the virtual source position to the iso-center 
and u is the distance from the central ray to the exposed detector column in fan 
direction. The line integrals are denoted by g and the derivative of the line integral with 
20 respect to a variable is expressed by the index. From the measured line integrals g 9 the 
line integrals g~ can be extrapolated for a virtual source position g according to 

g(u, a, a, q) = g(u, a, a,0) + Aqg~ (u, a, ar,0). (8) 

Therefore, Eq. (7) has to be solved for q . Reforming of Eq. (7) to 

= 2u (R 2 +u 2 ) J_ ua_ 

R 2 ^°^ Jf 2 £«a g&aa'*'~gr 8aa> (9) 

25 and partial integration with respect to u leads to 

(R 2 +u 2 ) 1 / ua V 

*» = " ^2 Sa +^ \gaa + — g aa ya, (10) 

which are the line integrals for the virtual source position q . 
The acquired and extrapolated data can now be resorted such that they 
correspond to an acquisition along a helical trajectory in the x-y-q space. Let R be the 
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vector from the center of rotation of the scanning system to the virtual source. The 
helical trajectory results in: 











R(a) = R 


sin or 


= R 


sin a 








< 4 J 



where a denotes the angular source position in relation to the x-axis. 
5 In a certain range, each value of q e [q^ , q m ] can be expressed by the 

linear equation 

q = l+T^l ^**™'* }^, (12) 

which fulfills the data acquisition of a helical trajectory in the x-y-q 
space. By this description, it is possible to define an offset a 0 as a starting point for the 

10 helical data acquisition, in order to use redundant data for the reconstruction process. 
This may lead to better image quality. 

The above described steps 1 to 4, in particular Step 3 may be applied and 
implemented in the CSCT scanner depicted and described with reference to Figs. 1 to 5, 
in the scanner depicted in Fig. 6 and in the data processing device depicted in Fig. 9. 

15 Fig. 8 shows an example for the characteristics of the wave-vector 

transfer q in dependence of the distance d of the scatter center from the detector. The 
distance S of the source to the center of rotation is chosen to be 570 mm, and the 
distance of the detector to the center of rotation is 470 mm, respectively. The distance R 
of the virtual source to the center or rotation is chosen to be 470 mm. The detector is 

20 supposed to consist of 20 rows of energy-resolving element with a size of 0.75 xO.75 
mm 2 in center of rotation. The x-ray source yields a spectrum of energies between 20 
and 160 keV. The object of 160 mm in diameter is placed in the center of rotation. The 
range of allowed ^-values is defined by the object size, the detector dimension and the 
energy spectrum. In the given example, q e [0.03,1.61], and all q values can be 

25 described by Eq. (12). The straight lines for the virtual source position at q = 1 .0 nm' 1 
are obtained by extrapolation according to Eq. 8. 

As may be taken from Fig. 8, characteristics of the wave-vector transfer 
in dependence of the distance d between the scatter center and the detector are depicted. 
For a given energy E, a certain detector element with a distance a to the scan plane or 
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slice plane receives radiation from the scatter centers arranged along a ray through the 
object under values being described by the hyperbolas. The straight lines are 
approximations to the corresponding hyperbola. In the range between and 
tfmax determined by E, h and the object diameter, every value for the wave-vector 

5 transfer may be expressed by equation (12). The straight lines intersect at Ihe virtual 
source position R. 

Fig. 9 shows an exemplary embodiment of a data processing device for 
performing steps 1 to 4 described above. As may be taken from Fig. 9, a central 
processing unit (CPU) or image processor 1 is connected to a memory 2 for storing 

10 read-outs from the detectors or the finally reconstructed data. As indicated before, the 
data may be acquired by a CSCT scanner such as depicted in Figs. 1 and 6. The image 
processor 1 may furthermore be connected to a plurality of input/output-network or 
other diagnosis devices. The image processor 1 is furthermore connected to a display 4 
(for example, to a computer monitor) for displaying information or images computed or 

1 5 adapted in the image processor 1 . An operator may interact with the data processor 1 via 
a keyboard 5 and/or other input or output devices, which are not depicted in Fig. 1 . 

The present invention described above, may, for example, be applied in 
the field of medical imaging. However, as described above, the present invention may 
also be applied in the field of non-destructive testing or baggage inspection. 
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